Orbital Kondo Effect in Cr0 2 : A LSDA+DMFT Study 
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Motivated by a collection of experimental results indicating the strongly correlated nature of the 
ferromagnetic metallic state of CVO2, we present results based on a combination of the actual 
bandstructure || with dynamical mean-field theory (DMFT) for the multi-orbital case. In striking 
contrast with LSDA(+U) J3| and model many-body approaches [j^), much better semiquantita- 
tive agreement with (i) recent photoemission results, (ii) domain of applicability of the half-metal 
concept, and (Hi) thermodynamic and dc transport data, is obtained within a single picture. Our 
approach has broad applications for the detailed first principles investigation of other transition 
metal oxide-based half-metallic ferromagnets. 



PACS numbers: 75.30.Mb, 74.80.-g, 71.55.Jv 

A lot of attention has recently been focused on CVO2, 
widely used in magnetic recording, and with potential ap- 
plications for spintronics |Q,0. In contrast to the CMR 
manganites, stoichiometric CrOi is already a ferromag- 
netic metal. Given the formal 4+ valence state of Cr, 
the two 3d electrons occupy t2 g orbitals. One would in- 
tuitively expect to form S — 1 spin on each site, and 
an antiferromagnetic Mott insulator. Why CVO2 is a 
ferromagnetic metal instead, has been answered by Ko- 
rotin et al. |J, who have carried out insightful (LDA 
+ U) calculations for this material. Their main conclu- 
sions are: (i) The triple-degeneracy of the ti g d-orbitals 
is lifted by tilting and rotation of the Cr — O octahedra in 
the basic rutilc structure, resulting in three bands with 
xy, yz ± zx character in the solid. The O 2p band(s) act, 
at least partially, as hole reservoirs resulting in Cr being 
mixed- valent (like Mn in doped manganites), explaining 
metallicity via self-doping in the negative charge-transfer 
situation realized in CrC?2- (H) an almost dispersionless 
majority spin band of predominantly d character at about 
1 eV below Ep over a large region of the Brillouin zone. 
This corresponds to strongly localized xy orbitals com- 
pletely occupied by one majority spin electron. But the 
d states (shown in Fig. ^) of predominantly yz ± zx char- 
acter hybridize with the O 2p band and disperse, cross- 
ing Ep. The Hund's rule coupling between the localized 
d xy spin and the spin density of the band d yz + zx elec- 
trons polarizes the latter, giving a ferromagnetic state 
via the double-exchange (DE) mechanism. Thus, both 
the metallicity and ferromagnetism are correlated well 
with each other. 

Closer examination reveals that CVO2 is a strongly 
correlated metal, implying that many-body correlation 
effects beyond the LDA Q (or its variants) need to be 
considered. A number of experimental observations sup- 
port such a view: 

(1) Polarization dependent XAS measurements reveal 
substantial ligand orbital polarization. An exchange 
splitting energy of A ex _ sp ; ~ 3.2 eV was deduced ||, 
implying substantial correlation effects, while LSDA cal- 



culations yield A ex _ sp i ~ 1.8eF! 

(2) The resistivity has a characteristic correlated Fermi 
liquid (FL) form §: p(T) = p Q + AT 2 + BT 7 / 2 @, 
and, in fact, CrO-x is a "bad metal" at high-T, with 
p dc (T > Tf M = 390K) exceeding the Mott limit [§. 
The Woods-Saxon ratio A/7 2 is close to that expected 
for heavy-fermion metals ]g] , implying substantial corre- 
lation effects in the Cr <i-band. 

(3) Optical conductivity studies reveal a Drude part at 
low energies, followed by a broad bump around 0.8 eV 
and high-energy features centered around 3 eV. LDA+U 
predicts only the small Drude part correctly ||. In ad- 
dition, noticeable spectral weight transfer (SWT) from 
high to low energy is found as T is reduced ||. This 
SWT scales with magnetization, M(T), as in the CMR 
materials Jlfj|, showing clearly the correlated nature of 
the metallic state (the SWT is a dynamical many body 
correlation effect, and cannot be accessed by LSDA+U). 

(4) Recent measurements show that the integrated 
photoemission (PES) lineshape is characterized by a 
low-energy quasicoherent feature along with an incoher- 
ent [[Fy|lJ] broad feature at lower energies. This satel- 
lite feature observed in PES is a signal for the impor- 
tance of dynamical, many-body correlation effects be- 
yond LDA+U. 

(5) Finally, recent optical |[o) and tunnelling mea- 
surements [H show half-metallic character only close 
to (about ~ 0.5 eV around) Ep. Direct comparison 
with LSDA+U shows that these calculations would yield 
half-metallicity up to 1.5 eV.lt is not easy to cure this 
within LSDA+U. One could, of course, refer to LSDA re- 
sults, where the minority-spin band does have a threshold 
around 0.6 eV, but LSDA gives results directly in conflict 
with (l)-(4). 

The emergence of a correlated FL scale required to 
understand the above features is out of reach of LSDA 
or pure DE models, because these are observed well be- 
low T C FM , and are thus related to additional scattering 
mechanisms in a half-metallic situation. We argued pre- 
viously UM that the above effects could be understood 
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by invoking the important role of local, dynamical or- 
bital correlations in the t 2g sector. However |Q, a model 
(gaussian) density of states was used there, limiting di- 
rect comparison to experimental results. To do this, one 
has to extend a model-based approach to include real 
bandstructure features via the LDA+DMFT method. 
We choose LDA+DMFT because although approxima- 
tions like LSDA+U do generate the correct ordered, 
insulating state(s), they fail to describe the correlated 
paramagnetic states, or the Mott transition accompa- 
nied by dynamical SWT. As discussed in detail pa] , very 
good agreement with experiment is achievable within 
LDA+DMFT. In practice, this is a highly non-trivial 
task, and very computationally expensive. Here, we de- 
scribe Cr0 2 within this "first-principles bandstructure" 
for correlated metals. The results are found to be in 
very good semiquantitative agreement with experiments 
probing both the occupied parts of the one-particle DOS 
(p(uj)) as well as with thermodynamic and dc transport 
data cited above. 

An understanding of features mentioned above should 
go hand-in-hand with the basic electronic structure. Our 
starting point is the LSDA+U work of Korotin et al. || , 
which yields the one-particle bandstructure of Cr0 2 (no- 
tice that the Hartree-Fock shift is already incorporated 
in LSDA+U). In Fig. [|, we show the partial t 2g DOS for 
the half-metallic situation realized in Cr0 2 . In the real 
material, only the t 2g states of Cr hybridized with 0-2p 
states are important, and so we neglect the higher energy 
e g bands. Following ||, we suppose the (3-2p states to 
introduce a small number of holes into the d yz ± zx bands. 
In fact, f*f 5 p\ 2s {uj)duj = 1.82 < 2. Given that the O- 

2p spectral density is small compared to p^ 2g (w), we keep 
only the t 2g bands in what follows. With the LSD A band- 
structure, the one-particle part of the Hamiltonian is: 



*/3( k ) c L c 



Ck/3o 



(1) 



k,a,/3 



with p T 2s (ui) = Y,k a p S(u) - e a/3 (k)) and a,(3 = xy, yz ± 
zx. To avoid double counting of Coulomb interaction 
contributions already contained in Hlda, a term H^ DA 
is subtracted from Hq. Following pa], the final result is, 



Ho=J2 e «/3( k ) c L<T C k/3<r + 5Z e ^ 

k,a,/3 iacr 



(2) 



where t\ aG = e iaa - U(n aa - 1/2) + (1/2) Jff(n a<7 - 1), 
and U and Jh are defined below in H . 

The full many body Hamiltonian for Cr0 2 is: 



H — H + - U^p n iaa nip a i — - JffS Jct -S 



i/3 



(3) 



To make progress, we notice that LSDA+U pushes the 
minority spin t 2g bands to high energy u > 2eV. To 



begin with, we focus on the majority-spin t 2g sector. A 
valid choice of parameters to describe this situation is 
given by the intra-orbital U = 5 eV, and the Hund's 
rule coupling Jh = 1 eV, close to the LSDA deduced 
values, along with the inter-orbital local interaction, U' = 
U — 2Jh, for t 2g orbitals. Notice that U', Jh couple the 
three xy,yz ± zx bands, requiring an extension to treat 
multi-orbital effects. 



1.2 



d. 



0.4 



0.0 




0.4 - 



3^ 
a.' 




1.2 



-2.0 



-1.0 



0.0 



1.0 



2.0 



3.0 



4.0 



CO 



FIG. 1. tig LSDA+U density of states per formula unit for 
both spin channels obtained from [3] , Figs. 3 (spin-down) 
and 4 (spin- up). 

We solve this multiband hamiltonian in d — 00 using 
the iterated perturbation theory (IPT), suitably gener- 
alized for the multi-orbital case. Notice that Cr0 2 falls 
into the class of TM oxides with well-separated t 2g and 
e g bands [[|. Moreover, G a p aa > (w) = 5 a p5 aa >G acr {uj) and 

^apae'i^) = 8 a p5 aa iY, at7 (uj). 

In the t 2g sub-basis, a DMFT solution involves (i) re- 
placing the lattice model by a self-consistently embedded 
multi-orbital, asymmetric Anderson impurity model, (ii) 
and a self-consistency condition which requires the local 
impurity Green function to be equal to the local GF for 
the lattice, given by 
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Vb 



d 6 k 
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(cj + p,)l-H° LSDA (k)-E(u;) 



(4) 



Using the locality of in d = 00, we have G a (uj) = 
G° a [uj — S Q (w)) using the Hilbert transform of the LSDA 
DOS. In contrast to |la], the inter-orbital coupling scat- 
ters electrons between the xy,yz ± zx bands, so only 
the total number, n^ t = Y]^ nd t2 , Q is conserved in a 
manner consistent with Luttinger's theorem. Further, in 
Cr0 2 , due to a highly asymmetric LDA DOS (Fig. [l]), 
we expect the final spectral function to reflect the inter- 
play of the complicated bandstructure and inter-orbital 
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correlations. 

The calculation follows the philosophy of the IPT for 
the one-band Hubbard model, but the self-energies and 
propagators are matrices in the orbital indices |l6[ . Leav- 
ing details for a longer paper, we present here the final 
equations. The local propagators are given by, 



AT 



1 



-/V ^— ' lu — X Q (u;) 

k 



(5) 



where a = xy,yz ± zx. The local self-energies are 
computed from a generalized IPT formalism that takes 
into account the Luttinger theorem constraint off half- 



filling 1 15 1 (generalized Friedel sum rule). Explicitly, 
A a [Ygl{u) + Sg(w) + E$M] 



E Q (w) = 



1 - B a [si 2 2(w) + T^Uu) + £$M] 



(6) 



where, for example, 7) 



E 

n.m 



(7) 



and G°(w) = [w + [i a - A^cj)]" 1 In Eqn. (| 

(l-27i Q )j7 et/j +e„-^ 



A = Mi^H and o 



a»(l-nj) 

Here, n a and n° are particle numbers determined from 
G a and G° a respectively, and D a ^[n] — (n aa nf3 a i) 
is calculated using (yi a o n po') = (n acr )(Ti/3cr') — 
Wtt I-oc /( a; )[Ea(w)G Q (a;)]c?a;. The last identity follows 
directly from the equations of motion for G a (cu). 
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t 2g DOS) of Cr0 2 for T = 0. Notice the high-energy 
"Hubbard bands" appearing as manifestations of dynam- 
ical effects of strong orbital correlations in the t 2g sector. 
As is also clear, the Luttinger theorem is obeyed to a 
very good accuracy. Additionally, the self-energies show 
characteristic correlated Fermi liquid behavior. Based on 
these results, we interpret the low energy quasicoherent 
features as arising from collective orbital Kondo screen- 
ing of tig orbital moments in the fully spin polarized 
half-metal. 

We emphasize the quantitative differences between the 
results of this work with an earlier one H], where an 
idealized (gaussian) unperturbed DOS was used; no low- 
energy (LSDA related) pseudogap-like feature was visible 
there. Further, the d xy band was replaced there by a dis- 
persionless, local level playing the role of polarizing the 
yz±zx bands via strong Jh, making it impossible to ac- 
cess changes resulting from the interplay of the realistic 
bandstructure and strong orbital correlations. These dif- 
ferences are especially important when one attempts to 

describe spectroscopy results quantitatively (see below). 
iui) 
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FIG. 3. The integrated photoemission lineshapes for the 
spin-up channel within LSDA+DMFT for U' = 3.0 eV (solid 
line) and LSDA+U (dashed line). 



0.3 



3 



0.6 



0.9 



-3.0 



-2.0 



-1.0 0.0 1.0 
CO 



2.0 



3.0 



4.0 



FIG. 2. LSDA+DMFT partial and total density of states 
for the Cr t2 g orbitals for U' — 3.0 eV. 

We now present our results. In Fig. |[ we show the 
renormalized, total spectral function (sum of the partial 



In Fig. ||, we show the integrated PES lineshape for low 
T, obtained directly from Ipes 

Quite satisfyingly, the quasicoherent spectral weight at \x 
is larger than what was reported in earlier experimental 
work |ll|], and in quite good agreement with results ob- 
tained recently from thin films In particular, it is 
gratifying to see good quantitative agreement between 
theory and experiment in the region —1.2 < uj < Ep. 
Specifically, the small, but clear dip seen in the LSDA+U 
DOS is absent in LSDA+DMFT, in agreement with the 
thin film results. At higher binding energy, we do not ex- 
pect quantitative agreement because both 0-2p and e g 
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bands of Cr begin to contribute, and also because the 
minority spin ti g band leads to a breakdown of our ap- 
proach, which neglected the minority-spin sector. Conse- 
quently, further work is needed to make comparison with 
the recent tunnelling and optical measurements |0,ff3l, 
since these measurements probe the minority spin band 
at higher energies. Optical measurements reveal that the 
almost complete spin polarization near Ep is reduced at 
higher energy u> > 0.5 eV. 

To this end, we extend the above calculation to in- 
clude the minority-spin t2 g bands, which (see Fig. are 
pushed to high energy oj > 2 eV by LSDA+U. Inter- 
orbital correlations between the tig T,| bands are ac- 
counted for in a way similar to that described above 
for the majority-spin sector. However, a larger XJ\i = 
V + 2 Jh is now required for the calculation of the mi- 
nority spin DOS, p^ 2fl (w), alongwith a Hartree shift of 
-Jhtti (where m is the magnetization per site) for the 
J,-spin sector. The interaction-corrected minority spin 
DOS (lower panel of Fig. ||) shows interesting features. 
First, U-[i broadens out the sharp LSDA features, and 
more importantly, transfers dynamical spectral weight 
from high- to lower energy u) ~ 0.7 eV . As a result, above 
0.7 eV, the spin polarization decreases continuously from 
saturation, in nice agreement with indications from op- 
tical and tunnelling measurements ^pUf, Notice that 
LSDA+U would predict complete half-metallicity up to 
a much higher energy ~ 1.5 eV. This represents strong 
evidence for the importance of treating dynamical SWT 
correctly, an effect missed by LSDA+U, but treated ad- 
equately by LSDA+DMFT. In fact, in our view, opti- 
cal and tunnelling results provide a direct confirmation 
of the importance of dynamical correlation effects. In- 
deed, one could argue that the near similarity of the 
LSDA- and LSDA+DMFT spectral functions for t\ g sec- 
tor for —0.7 < oj < 0.2 eV is a sign of the adequacy 
of LSDA+U. As discussed above, however, consistency 
with optics and tunnelling results as well is only achieved 
within LSDA+DMFT. 

The results obtained above are also consistent with the 
low-T thermodynamic and dc resistivity of CrOi- In- 
deed, from the computed self-energy, we estimate a mod- 
erate quasiparticle renormalization, Z. From the fact 
that ImSy Z ± zx (oj) ~ —bto 2 , we infer that the low-T dc re- 
sistivity should follow pdc{T) = pa + AT 2 . In fact, using 
A = (m* /ne 2 ){d 2 Y.{uj)/duj 2 ) u= E F from Im^ yz±zx (uj), 
we estimate the Woods-Saxon ratio, A/-f 2 ~ 3.5 10~ 5 , 
close to the value found experimentally ||] . 

To conclude, we have extended the LSDA+U calcu- 
lation || to explicitly include the dynamical correla- 
tion effects arising from local electronic correlations in 
the tig sector in the multi-orbital system CrOi- Puz- 
zling signatures of strong correlations in the half-metallic 
state are understood as manifestations of collective or- 
bital Kondo effect in the ti g sector, originating from a 
non-trivial interplay between the realistic hopping in the 



rutile structure and inter-orbital (t2 g ) correlations. Much 
better quantitative agreement with PES is obtained with 
LSDA+DMFT than with LSDA or LSDA+U, showing 
clearly the importance of dynamical correlation effects. 
Further, our results are also in semiquantitative agree- 
ment with thermodynamic and dc transport measure- 
ments. Essentially similar techniques can also be used 
fruitfully for other transition metal oxide-based ferro- 
magnets currently of great interest (l^] . 

Work carried out under the auspices of the Sonder- 
forschungsbereich 608 of the Deutsche Forschungsge- 
meinschaft. 
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